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We present a novel solution method with which we obtain the spectral functions of the attractive Hubbard model 
in two dimensions in a non self-consistent formulation of the T-matrix approximation. We use a partial fraction 
decomposition of the self energy from which we obtain the single-particle spectral functions to a remarkably high 
accuracy. Here we find the poles and residues of the self energy to a relative accuracy of 10 -80 for an 8 x 8 lattice, 
and plot the resulting density of states. Our results show pseudogap physics as T c is approached from above; 
these analytical results are in agreement with recent Quantum Monte Carlo data. 



One can motivate a study of the attractive 
Hubbard model (AHM) by suggesting that it 
might provide a rudimentary description of the 
correlated electrons in the high-T c materials. 
One feature that theorists are keen to see re- 
produced in any model is the appearance of a 
normal-state pseudogap |jj above T c that evolves 
smoothly into the energy gap of the supercon- 
ducting state. For example, QMC studies |2) have 
produced evidence of pseudogap physics in the 
two-dimensional AHM. 

If the T-matrix approximation applied to the 
AHM can be shown to suppress the density of 
states (DOS) at the Fermi level for temperatures 
above T c , this would lend credence to the notion 
that the pseudogap is the result of strong pair- 
ing fluctuations. In this brief report, we demon- 
strate that one can calculate the DOS to an arbi- 
trarily high accuracy for the non self-consistent 
(NSC) version of the T-matrix approximation, 
and that one in fact obtains results very simi- 
lar to the QMC data. We stress that this calcu- 
lational method should be applicable to any non 
self-consistent microscopic theory involving mero- 
morphic functions. 

To begin, we must calculate the fully interact- 
ing DOS within the NSC T-matrix approxima- 
tion. Such a calculation is made possible by the 
fact that we have access to the exact solution of 
the NSC pair propagator, viz. 



G^_(Q^n) = X(Q>«)/(1 - \U\x(Q,iVn)) 

where the analytic continuation of x{Qi^ v n) to 
the complex plane is given by 
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The latter function is meromorphic with a finite 
number of simple poles on the real axis (as usual, 
/ is the Fermi-Dirac distribution and £ is the non- 
interacting single-particle energy relative to the 
chemical potential). Likewise, the analytic con- 
tinuation of the pair propagator, 
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can be shown to be a meromorphic function with 
a finite number of simple poles. This demon- 
strates that the pair propagator can be written 
as a rational polynomial with real coefficients, a 
consequence of which is that there are only two 
possibilities for the placement of its singularities: 
they lie either on the real axis (corresponding to 
the normal state) or in conjugate pairs equally 
spaced above and below the real axis (correspond- 
ing to the unstable state below T c ). Moreover, 
since the degree of the polynomial in the denomi- 
nator is larger by one than that of the polynomial 
in the numerator, the high-frequency asymptotic 
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behaviour of the pair propagator is governed by 
G v +-{Q,z) ~ 1/z as \z\ -> oo. 

Thus, in the normal state, the pair propaga- 
tor admits a partial fraction decomposition into 
a series of simple poles. We write this as 
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where the energies £V are real and the residues 

> are strictly positive. [Below we denote 

the number of such poles for each Q component 
by Sq-] Quite simply, such a decomposition can 
be executed numerically in MapleV to arbitrary 
accuracy (without an associated increase in com- 
puting time), and in what follows we used Dig- 
its:=80 to obtain a relative accuracy of 10 -80 . 

We now show how this result can be exploited 
to give us a practical, numerical method for cal- 
culating the self-energy. The first step is to put 
the partial fraction form of the pair propaga- 
tor, Eq. (0), into the self energy using the iden- 
tity G^_(Q, iv n ) = G^_(Q,z = iv n ) to recover 
the Matsubara frequency components. Then, we 
complete various Matsubara frequency sums, and 
the final result is 
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where b is the Bose distribution function. Once 
the self-energy is calculated numerically by this 
procedure, the Green's function follows from 
Dyson's equation, and the single-particle spectral 
function follows in the usual way. The DOS can 
then be constructed according to 
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This yields a DOS that is a series of ^-function 
peaks; when plotting our results we incorporate a 
small artificial broadening to smooth the curves. 
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Above we show two figures for the resulting 
DOS for an 8 x 8 two-dimensional square lattice 
for \U\ being half the band width. The chemical 
potential is varied to fix the density to be that 
of a 1/4- filled system as the temperature is low- 
ered. The pseudogap is clearly visible, whereas 
at T/t = 1 (not shown) the spectrum shows no 
anomaly at the Fermi level. Note that in this 
system (fixed density ^ one) T c = ||. 
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